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HIGHLIGHTS 


•  Pyrolysis  of  coal,  straw,  and  blended  coal-straw  particles  were  studied  by  using  a  single-particle  reactor  system. 

•  Experimental  data  indicated  a  lack  of  synergistic  effects. 

•  A  model  was  developed  based  on  the  CPD  model  to  simulate  the  pyrolysis  process  of  coal  and  biomass  particles. 

•  An  extended  CPD  model  was  proposed  to  describe  the  co-pyrolysis  characteristics  of  coal-biomass  blends. 
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This  paper  reports  a  combined  experimental  and  numerical  investigation  of  the  pyrolysis  characteristics 
of  coal,  biomass,  and  coal-biomass  blends.  Coal  and  straw  were  grounded  and  pressed  into  spherical  par¬ 
ticles  with  diameter  of  8  mm,  and  blended  coal-straw  particles  were  prepared  through  mixing  pulverized 
coal  and  straw  before  pressed  into  particles.  Sample  particles  were  suspended  in  the  center  of  a  single¬ 
particle  reactor  system  and  devolatilized  under  different  temperatures.  The  analysis  of  the  time  history  of 
the  residual  mass  of  particles  of  coal,  straw,  and  coal-straw  blends  suggested  an  absence  of  synergistic 
effect  between  the  coal  and  the  straw.  In  addition,  a  one-dimensional,  time-dependent  particle  model; 
based  on  the  chemical  percolation  devolatilization  (CPD)  and  bio-CPD  models;  was  developed  to  simulate 
the  pyrolysis  of  coal  and  straw  particles.  The  model  predictions  agree  will  with  the  measured  data.  An 
extended  CPD  model  was  proposed  to  explain  the  co-pyrolysis  characteristics  of  coal-biomass  blends. 
Encouraging  agreement  was  found  between  the  predicted  and  the  experimental  results  of  pyrolysis  of 
coal-straw  blends. 

©  2014  Elsevier  Ltd.  All  rights  reserved. 


1.  Introduction 

With  the  growing  shortage  of  fossil  fuels  around  the  world, 
development  and  utilization  of  biomass  energy  has  attracted  lots 
of  attentions.  As  a  sustainable  fuel,  biomass  generates  energy  from 
biological  material  like  straw,  sawdust,  wood  waste,  etc.  Co-firing 
of  coal  and  biomass  fuels  is  considered  as  an  alternative  effective 
method  for  coal-fired  power  plants,  since  it  helps  to  reduce  the 
emission  of  C02  and  other  greenhouse  gases  [1-3  .  However,  the 
application  of  this  technology  in  power  plant  boilers  has  many 
technical  challenges  [4-6].  Pyrolysis  can  convert  low  chemical 
energy  density  fuels,  such  as  biomass  and  low  rank  coals,  to  high 
energy  density  fuels,  such  as  gas,  and  char.  But  the  pyrolysis 
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kinetics  of  coal-biomass  blends  is  still  not  clear,  which  is  extre¬ 
mely  important  to  understand  the  co-combustion  processes. 

Co-pyrolysis  of  coal  and  biomass  mixtures  has  been  investi¬ 
gated  by  several  researchers.  Most  previous  studies  [7-11  support 
that  there  is  no  synergistic  effect  between  coal  and  biomass  during 
the  co-pyrolysis  process,  except  that  a  few  studies  [12]  find  the 
existence  of  a  synergy.  Thermo-gravimetric  analysis  (TGA),  which 
measures  the  mass  loss  of  samples  as  a  function  of  temperature, 
is  the  most  common  technique  utilized  to  study  pyrolysis  of 
coal-biomass  mixtures  [7-10].  However,  in  these  studies  the  heat¬ 
ing  rate  often  varies  between  10  K/min  and  50  K/min,  due  to  the 
limitation  of  TGA.  It  has  been  reported  that  the  heating  rate  has 
an  important  effect  on  the  pyrolysis  reaction  kinetics  and  products 
yield  [13-15].  The  kinetic  parameters  derived  from  pyrolysis 
experiments  at  low  and  high  heating  rates  are  quite  different 
[13].  Many  researches  [15-17]  have  shown  evidence  of  increase 
in  the  volatiles  yield  as  the  heating  rate  elevated.  In  the  experiment 
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Nomenclature 

Cp 

specific  heat,  J / ( kg  I<) 

Cl 

first  radiation  constant 

Subscrivts 

C2 

second  radiation  constant 

0 

initial  condition 

d 

diameter,  m 

P 

particle 

Ebz 

black  body  emissive  power 

qz 

quartz 

h 

convective  heat  transfer  coefficient,  W/(m2  K) 

s 

solid 

Nu 

Nusselt  number 

V 

volatile 

Pr 

Prandtl  number 

w 

furnace  wall 

r 

radius  coordinate,  m 

oo 

nitrogen  in  the  reactor 

RP 

particle  outer  radius,  m 

2 

wavelength 

Re 

Reynolds  number 

t 

time,  s 

Definition  of  the  kinetic  parameters  for  the  CPD  model 

T 

temperature,  K 

Eb 

bridge  breaking  activation  energy,  kcal/mol 

u 

velocity  of  volatile,  m/s 

bridge  breaking  frequency  factor,  s  1 

Cb 

standard  deviation  of  Eb,  kcal/mol 

Greek  symbols 

Eg 

gas  formation  activation  energy,  kcal/mol 

P 

density,  kg/m 

Ag 

gas  release  frequency  factor,  s_1 

A 

thermal  conductivity,  W/(m  I<) 

°g 

standard  deviation  of  Eg,  kcal/mol 

CO 

volatile  generation  rate,  s_1 

P' 

kinetic  ratio  of  bridge  breaking  to  char  formation 

AH 

heat  of  pyrolysis  reaction,  kj/kg 

Ec 

difference  in  activation  energy  between  bridge  breaking 

a 

Stefan-Boltzmann  constant,  W/(m2  K4) 

and  char  formation,  kcal/mol 

s 

emissivity 

F 

Across 

cross-linking  activation  energy,  kcal/mol 

S 

thickness,  m 

kcross 

cross-linking  frequency  factor,  s_1 

a 

absorptivity 

z 

transmissivity 

y 

reflectivity 

by  Gibbins-Matham  and  Kandiyoti  [17],  the  yield  of  total  volatiles 
on  a  dry  ash-free  (daf)  basis  increased  from  47%  at  1  K/s  to  52%  at 
1000  K/s  with  coal  samples  heated  to  a  final  temperature  of  973  K 
and  held  for  30  s.  In  this  work,  a  single-particle  reactor  system  was 
used  to  measure  the  time  history  of  the  mass  loss  and  temperature 
rise  of  samples  simultaneously,  with  a  significantly  higher  heating 
rate  (more  than  1000  K/min)  than  TGA. 

Numerous  kinetic  models  with  various  reaction  kinetic  param¬ 
eters  for  pyrolysis  of  coal  and  biomass  have  been  reported  [18-20]. 
These  different  kinetic  models  play  an  important  role  in  previous 
particle  modeling  researches  [20-22].  Among  these  models,  the 
chemical  percolation  devolatilization  (CPD)  model  [23]  is  distinc¬ 
tive,  since  its  kinetic  rates  are  general,  not  specific  to  a  certain  type 
of  coal  or  biomass  and  valid  over  a  wide  range  of  temperatures  and 
heating  rates  [24].  In  previous  studies,  the  CPD  and  bio-CPD  (an 
upgraded  CPD)  models  can  accurately  predict  the  pyrolysis  process 
of  different  kinds  of  pulverized  coal  and  biomass  in  different  reac¬ 
tors  [23,25-27  .  Despite  these  efforts,  studies  on  a  general  pyroly¬ 
sis  model  of  the  pyrolysis  of  blended  coal-biomass  particles  have 
been  rarely  reported.  In  this  study,  the  CPD  and  bio-CPD  models 
were  developed  further  to  describe  the  pyrolysis  of  large  coal 
and  biomass  particles;  and  the  co-pyrolysis  of  blended  coal-bio¬ 
mass  particles  was  modeled  by  using  an  extended  CPD  model. 

The  objective  of  this  study  was  to  investigate  the  co-pyrolysis 
characteristics  of  coal-biomass  blends  under  high  heating  rates 
(compared  to  TGA),  by  recording  the  time  history  of  the  mass  loss 
and  temperature  rise  of  a  single  particle;  and  to  develop  a  general 
pyrolysis  model  of  coal-biomass  blends. 

2.  Experimental  setup 

Pyrolysis  experiments  were  conducted  for  a  single  coal,  straw 
and  blended  coal-straw  particle.  Coal  and  straw  were  grounded 
and  pressed  into  spherical  particles  with  diameter  of  8  mm,  and 
blended  coal-straw  particles  were  made  through  mixing 


pulverized  coal  and  straw  in  different  proportions  before  the  press¬ 
ing  process.  All  particles  were  dried  at  378  I<  for  30  min  before  the 
pyrolysis  experiments  to  eliminate  the  influence  of  water  vapori¬ 
zation.  The  coal  sample  used  is  a  Chinese  brown  coal  called  Zhun- 
dong  brown  coal.  Its  proximate  analysis,  ultimate  analysis  and 
biomass  characterization  are  listed  in  fable  1. 

Pyrolysis  experiments  were  carried  out  in  a  single-particle  reac¬ 
tor  system,  which  is  shown  in  Fig.  1.  The  system  was  composed  of 
three  main  parts:  the  quartz  reactor,  the  heating  furnace,  and  the 
mass  and  temperature  data  loggers.  A  single  coal,  biomass  or 
blended  coal-biomass  particle  was  suspended  on  type-S  thermo¬ 
couple,  which  was  inserted  into  the  center  of  the  particle.  Both 
the  thermocouple  and  the  particle  were  rested  on  a  mass  sensor 
to  simultaneously  obtain  the  mass  and  temperature  at  the  particle 
center.  A  data  acquisition  system  (Model  Agilent  34970A)  recorded 
temperature  data  from  the  thermocouple  at  1  Hz.  The  mass  sensor 
(Model  Beijing  Hengjiu  Instrument  HTG-1)  recorded  the  time  his¬ 
tory  of  particle  mass  loss  at  a  resolution  of  0.01  mg  and  an  accuracy 
of  ±0.1  mg. 


Table  1 

Chemical  analyses  of  the  coal  and  straw  samples. 

Ash  Volatile  Fixed  carbon 

Proximate  analysis  (dry  basis  wt%) 

Coal  4.34  30.86  64.79 

Straw  10.08  72.04  17.87 


C  H  N  S  O 

Ultimate  analysis  (dry  basis  wt%) 

Coal  75.39  3.48  1.19  0.42  15.19 

Straw  46.65  6.01  1.01  0.20  36.05 

Cellulose  Hemicellulose  Lignin 

Biomass  characterization  (dry  basis  wt%) 

Straw  39.89  37.77  21.94 
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Fig.  1.  Schematic  diagram  of  the  single-particle  reactor  system. 


An  experiment  began  by  suspending  the  thermocouple  and  par¬ 
ticle  on  the  mass  sensor.  Then  the  quartz  reactor  was  installed 
through  the  bottom  of  the  thermocouple.  The  electric  furnace 
was  heated  to  preset  temperature  and  then  lifted  up  quickly  to 
the  height  of  the  quartz  reactor.  The  particle  was  heated  by  radia¬ 
tion  of  the  furnace  wall  and  its  temperature  and  mass  loss  data 
were  recorded  simultaneously  once  every  one  second.  All  experi¬ 
ments  were  conducted  in  nitrogen  with  the  flow  rate  of  5  L/min 
(20  °C,  1  atm),  in  which  particles  experienced  the  pyrolysis 
process. 

The  advantage  of  this  experimental  setup,  as  compared  to  con¬ 
tinuous  feeding  setups,  e.g.  [26],  is  the  ability  to  obtain  the  tem¬ 
perature  and  mass  loss  data  of  a  particle  during  the  pyrolysis 
process  simultaneously,  which  is  extremely  important  for  the 
understanding  of  the  pyrolysis  characteristics  and  the  validation 
of  models. 


3.  Pyrolysis  model 


precursors  are  cross-linked  back  into  the  char  matrix,  which  can¬ 
not  vaporize  because  the  molecular  weights  are  too  large.  The  def¬ 
initions  of  the  CPD  kinetic  parameters  are  shown  in  the 
nomenclature  table. 

The  chemical  structural  parameters  used  in  the  CPD  model  are 
the  molecular  weight  of  the  cluster  (MWd),  the  molecular  weight  of 
side  chains  ( MWS ),  the  initial  fraction  of  intact  bridges  (p0),  the 
coordination  number  (cr  +  1 ),  and  the  initial  fraction  of  char  bridges 
(c0).  The  first  four  parameters  are  typically  determined  from  13C 
NMR  measurements,  while  the  fifth  c0  cannot  be  measured  directly 
and  was  determined  empirically. 

The  bio-CPD  model,  developed  by  Fletcher  et  al.  [25]  and  other 
researchers  [26,27,29],  is  used  to  model  biomass  pyrolysis  in  this 
work.  Considering  that  the  CPD  model  can  be  used  for  components 
of  biomass  [25  ,  which  is  similar  to  a  low  rank  coal,  the  bio-CPD 
model  calculates  the  pyrolysis  of  biomass  as  a  linear  combination 
of  the  constituent  compounds,  i.e.,  cellulose,  hemicellulose,  and 
lignin.  The  amounts  of  these  compounds  are  normalized  so  that 
they  sum  to  100%  [26  .  Based  on  the  bio-CPD  model,  this  work 
extends  the  CPD  model  to  model  the  co-pyrolysis  of  coal-biomass 
blends  by  combining  the  predicted  pyrolysis  yields  of  the  coal  and 
biomass.  The  chemical  structural  parameters  and  kinetic  parame¬ 
ters  of  the  model  are  summarized  in  Tables  2  and  3  from  the  liter¬ 
atures,  except  that  the  structural  parameters  of  the  coal  are 
estimated  by  using  the  nonlinear  modified  quadratic  correlation 
of  13C  NMR  measurements  of  volatile  matter  content  and  coal 
structure  with  ultimate  analysis  [30  . 

The  CPD  model  assumes  that  the  particle  is  isothermal,  which  is 
only  acceptable  when  the  particle  is  sufficiently  small.  In  general, 
when  the  Biot  number  is  less  than  0.1  [31],  the  intraparticle  heat 
transfer  is  negligible  and  the  CPD  model  can  be  applied.  Flowever, 
to  model  the  pyrolysis  of  a  large  particle,  an  energy  equation  must 
be  used  to  predict  the  temperature  inside  the  particle.  Coupling  the 
CPD  model  and  energy  equation,  the  time  history  of  the  tempera¬ 
ture  and  mass  loss  during  the  pyrolysis  of  a  large  particle  can  be 
simulated. 

The  particle  energy  equation  is  based  on  the  one-dimensional 
unsteady  heat  conduction  equation  in  a  spherical  coordinate, 
which  assumes  that  the  particle  is  spherically  symmetric: 

dT  1  d  2  1  d  f  2  dT\ 

pCp*  Yt  +  Pv^v  ~p2  Y (r  Tu)  =  Wd?(rAsdi:)+coAH  (1 } 


The  nomenclature  table  provides  a  complete  listing  of  the  sym¬ 
bols  and  subscripts. 

Under  the  assumption  of  immediate  outflow  of  volatiles  and 
thermal  equilibrium  between  volatiles  and  solid  in  any  control  vol¬ 
ume  [22],  the  velocity  u  of  volatiles  passing  across  a  control  surface 
at  a  distance  r  from  the  particle  center  is  presented  in  Eq.  (2). 


Fletcher  et  al.  [23,28]  developed  the  chemical  percolation  dev¬ 
olatilization  (CPD)  model  to  predict  coal  pyrolysis  yields  as  a  func¬ 
tion  of  time,  using  a  description  of  the  unique  chemical  structure  of 
various  coals.  Coal  is  modeled  as  a  series  of  aromatic  clusters  con¬ 
nected  by  labile  bridges.  All  coals  use  the  same  kinetic  parameters. 
During  the  heating  of  coal  pyrolysis,  the  activating  and  breaking  of 
bridges  are  modeled  with  two  competing  pathways.  The  activated 
intermediate  bridges  can  either  break  into  side  chains  or  form  sta¬ 
ble  char  bridges.  The  two  competing  reaction  rates  are  a  function  of 
kinetic  parameters,  and  the  relationship  between  the  number  of 
broken  bridges  and  detached  clusters  is  predicted  using  percola¬ 
tion  statistics  for  Bethe  lattices.  The  tar  and  gas  yields  are  not  input 
parameters  but  calculated  using  the  kinetic  mechanism,  a  flash  cal¬ 
culation,  and  a  vapor  pressure  correlation.  In  addition,  parts  of  tar 


w  =  _d(pM) 

dt 


Table  2 

Structural  parameters  to  model  coal  and  biomass  pyrolysis  using  the  CPD  model. 


Structural  parameter 

Coal 

Cellulose 

Hemicellulose 

Lignin 

MWC/ 

329 

162 

162 

197 

MW,, 

25 

37 

37 

37 

Po 

0.70 

1.00 

1.00 

0.71 

a  +  1 

5.6 

3.0 

3.0 

3.5 

c0 

0.35 

0.00 

0.00 

0.00 

Ref. 

[26,29] 

[26,29] 

[26,29] 

I<.  Wan  et  al./Fuel  139(2015)  356-364 


359 


Table  3 

Kinetic  parameters  to  model  coal  and  biomass  pyrolysis  using  the  CPD  model  [25-27]. 


Kinetic  parameter 

Coal 

Cellulose 

Hemicellulose 

Lignin 

Et,  (kcal/mol) 

55.4 

54.1 

47.6 

54.0 

/Ms-') 

2.6  x  1015 

2.1  x  1015 

8.0  x  1014 

2.6  x  1015 

Gb  (kcal/mol) 

1.8 

2.7 

1.9 

4.0 

Eg  (kcal/mol) 

69.0 

61.2 

38.2 

69.0 

3.0  x  1015 

3.0  x  1015 

3.0  x  1015 

2.3  x  1019 

ag  (kcal/mol) 

8.1 

8.1 

5.0 

2.6 

P' 

0.9 

3.0 

1.6 

3.9 

Ec  (kcal/mol) 

0.0 

0.0 

0.0 

0.0 

Ecmss  (kcal/mol) 

65.0 

65.0 

65.0 

65.0 

Across  (S  ) 

3.0  x  1015 

3.0  x  1015 

3.0  x  1015 

3.0  x  1015 

Eq.  (3)  is  the  conservation  of  mass  for  the  particle.  The  volume 
of  the  solid  is  assumed  to  be  constant  while  its  density  decreases  in 
proportion  to  the  total  volatile  matter  loss  during  the  pyrolysis 
process  [20,21  .  The  pyrolysis  rate  co  is  calculated  by  the  CPD 
model,  which  obtains  the  temperature  of  each  control  volume  from 
the  energy  equation.  Thus,  the  temperature  and  mass  loss  of  parti¬ 
cles  can  be  predicted  simultaneously. 

The  surface  of  a  particle  is  subjected  to  both  connective  and 
radiative  heat  flux,  which  can  be  modeled  by: 

f)T 

-as  —  \r=Rp  =  hp(TRp  -  T„)  +  %£p<r(7tp  -  K)  +  Zp°(T4rp  ~  O 

(4) 

The  convective  heat  transfer  coefficient  of  the  particle,  hpi  is 
based  on  the  correlation  [21]: 

=  2  +  0.6Re1/2Pr1/3  (5) 

■^oo 

The  spectral  transmission  rate  of  quartz  reactor  iqZtA  and  black 
body  emissive  power  Ehx  under  different  furnace  wall  tempera¬ 
tures  is  shown  in  Fig.  2.  It  can  be  found  that  the  quartz  reactor 
transmits  radiation  partly  in  shorter  wavelength  region  than 
3.7  pm.  The  transmissivity  of  quartz  reactor  Tqz  is  presented  in 
Eq.  (6). 

fo  Tqz^EbldA  C\X* 

Tqz  bk  ec2/(A*Tw)  _  1 

u  1  w 


Temperature  of  quartz  wall  can  be  modeled  by  using  Eq.  (7),  the 
heat  transfer  equation  of  quartz  reactor  wall.  The  radiation  absorp¬ 
tivity  of  quartz  reactor  is  defined  by  Eq.  (8)  and  the  convective  heat 
transfer  coefficient  of  quartz  reactor  is  presented  in  Eq.  (9). 

cfT 

-CP,qzPqzSv^f  =  MT, z  -  Too)  +  0^0(1%  -  T4J  (7) 


0 

O 

CL 

0 

> 

CO 

CO 

E 

LU 


Fig.  2.  Spectral  transmission  rate  of  the  quartz  reactor  and  black  body  emissive 
power  under  different  furnace  wall  temperatures. 


NuqzA00 

Dqz 


The  model  requires  a  full  set  of  physical  properties  of  the  parti¬ 
cle  and  the  quartz  reactor,  which  is  summarized  in  Table  4  and  is 
based  on  practical  values  or  recommendations  from  literatures. 
Most  physical  properties  of  the  biomass  are  set  to  be  the  same  as 
those  of  the  coal  as  a  reasonable  simplification,  except  for  the  pyro¬ 
lysis  heat. 


4.  Results  and  discussion 

In  this  section,  pyrolysis  data  (particle  temperature  vs.  time  and 
residual  mass  vs.  time)  of  coal,  straw,  and  blended  coal-straw  par¬ 
ticles  were  collected  using  the  single-particle  reactor  system  and 
compared  with  model  predictions.  With  these  data,  the  pyrolysis 
characteristics  of  coal-straw  blends  were  systematically  analyzed. 

4.1.  Pyrolysis  of  coal  particles 

Figs.  3  and  4  illustrate  the  time  history  of  the  temperature  and 
residual  mass  of  coal  particles  under  different  furnace  wall  temper¬ 
atures  (1000  K,  1100  K,  and  1200  K),  respectively.  All  particles  had 
the  similar  initial  volume  (diameter  of  8  mm)  and  initial  mass 
(~360  mg).  Each  experiment  was  repeated  three  times  and  the 
average  result  is  shown  together  with  the  error  bar  indicating 
the  statistical  uncertainty  of  the  three  measurements  at  each  con¬ 
dition.  As  plotted  in  Fig.  3,  at  the  beginning  of  the  pyrolysis  the  coal 
particle  temperature  maintained  at  300  I<  (the  ambient  tempera¬ 
ture)  for  a  few  seconds,  since  the  heat  had  not  been  conducted  into 
the  center  of  the  particle.  Then  the  temperature  rose  rapidly  due  to 
the  radiative  heating  of  the  furnace  wall  and  reached  a  stable  value 
in  the  final  stage  of  the  pyrolysis,  which  means  the  particle 
achieved  thermal  equilibrium  in  the  furnace.  From  Fig.  4,  it  can 
be  found  that  the  particle  gradually  lost  its  weight  as  the  temper¬ 
ature  rose  and  volatile  yielded  during  the  pyrolysis  process.  In  the 
final  stage,  the  mass  loss  of  the  coal  particle  slowed  down  and  the 
residual  mass  remained  stable  finally. 

At  the  time  of  500  s,  the  measured  total  volatile  yields  of  the 
coal  particle  were  25.9%,  30.9%,  and  31.6%  (dry  basis)  with  the  fur¬ 
nace  wall  temperature  at  1000  K,  1100  K,  and  1200  K,  respectively. 
With  the  wall  temperature  rising  from  1000  I<  to  1100  K,  the  vola¬ 
tile  yields  of  the  particle  increased  significantly  (5%),  whereas  the 
increase  of  volatile  yields  became  less  (0.7%)  with  the  wall  temper¬ 
ature  continually  rising  to  1200  K.  Since  the  volatile  mass  fraction 
from  the  proximate  analysis  of  the  coal  was  30.86%  (dry  basis; 
measured  at  1173  K;  see  ^able  1),  it  can  be  concluded  that  the  coal 
particle  almost  reached  the  maximum  yields  at  1100  I<  and  higher 
temperature  could  not  enhance  the  volatile  yields  significantly. 

By  comparing  the  experimental  data  with  model  predications,  it 
can  be  found  that  the  CPD  model  can  reasonably  predict  the  pyro¬ 
lysis  of  coal  particles  under  different  temperatures.  Except  that  the 
modeled  pyrolysis  at  1200  K  was  slightly  deviated,  the  prediction 
data  showed  excellent  agreement  with  the  experimental  data.  At 
the  time  of  500  s,  the  model  predicted  that  the  total  volatile  yields 
of  the  coal  particle  were  26.5%,  31.2%,  and  32.1%  (dry  basis), 
respectively,  at  the  three  furnace  wall  temperatures.  The  predicted 
volatile  yields  matched  the  measured  data  within  0.6  wt%  upon 
500  s  of  the  pyrolysis  process. 

Fig.  5  shows  the  time  history  of  the  mass  loss  rate  of  the  coal 
particle  under  different  furnace  wall  temperatures.  The  measured 
maximum  pyrolysis  rate  increased  from  0.19%  s-1  (1000  K)  to 
0.32%  s”1  (1100K)  and  finally  reached  0.49%  s”1  (1200  I<),  as  the 
furnace  wall  temperature  rose.  And  the  time  corresponding  to 
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Table  4 

Physical  properties  of  the  particle  and  the  quartz  reactor. 


Variable 

Value 

Ref. 

d 

Diameter  of  particle,  m 

8.0  x  10“3 

m0 

Initial  mass  of  particle,  kg 

360  x  10“6 

sp 

Emissivity  of  particle 

0.9 

[32] 

Cp,s 

Specific  heat  of  solid,  J / ( kg  K) 

1150 

for  T  <  573  K 

[33] 

1150  +  2.03(1-573)  -  1.55 

x  lO~3(T-  573)2  for  T  >  573  K 

Cp,V 

Specific  heat  of  volatile,  J / ( kg  K) 

1100 

[22] 

4 

Thermal  conductivity  of  particle,  W/(m  K) 

0.19 

for  T  <  573  K 

[34] 

0.19  +  2.5  x  10^4(T  -  573)2 

for  T  >  573  K 

AH 

Heat  of  pyrolysis  reaction,  kj/kg 

418.6  for  coal 

[35] 

210.0  for  biomass 

[36] 

Dqz 

equivalent  diameter  of  quartz  reactor,  m 

26.0  x  10“3 

dqz 

Thickness  of  quartz  reactor,  m 

2.0  x  10“3 

Cp.qz 

Specific  heat  of  quartz,  J /( kg  K) 

931.3  +  0.256T-24.0T~2 

[37] 

Pqz 

Density  of  quartz,  Kg/m3 

2200 

[38] 

y  qz 

Reflectivity  of  quartz 

0.08 

[38] 

Fig.  3.  Time  history  of  center  temperature  of  coal  particles  under  different  furnace 
wall  temperatures. 


Fig.  5.  Time  history  of  mass  loss  rate  of  coal  particles  under  different  furnace  wall 
temperatures. 


Time  (s) 

Fig.  4.  Time  history  of  residual  mass  of  coal  particles  under  different  furnace  wall 
temperatures. 

the  maximum  pyrolysis  rate  decreased  from  93  s  to  46  s.  Hence, 
the  higher  the  pyrolysis  temperature  coal  particles  experienced, 
the  larger  the  maximum  pyrolysis  rate  became  and  the  less  the 
time  is  needed.  Compared  with  the  measured  data,  the  CPD  model 
tended  to  overestimate  the  maximum  pyrolysis  rate  and  underes¬ 
timate  the  time  to  reach  the  maximum  rate.  However,  except  for 


Time  (s) 

Fig.  6.  Time  history  of  mass-weighted  average  temperature  (lines  with  symbols) 
and  temperature  rise  rate  (lines)  of  coal  particles  under  different  furnace  wall 
temperatures. 


the  disagreement  of  these  values,  the  model  predicted  the  trends 
of  the  mass  loss  rate  and  its  variation  with  temperature  well. 

The  mass-weighted  average  temperature  of  coal  particles  and 
the  rise  rate  of  the  temperature  calculated  by  the  CPD  model  are 
shown  in  Fig.  6.  The  mass-weighted  average  temperature  was  cal¬ 
culated  by  averaging  the  temperature  of  each  node  with  the  mass 
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of  the  node  as  the  weighting  factor,  and  can  be  used  as  a  character¬ 
istic  temperature  for  the  whole  particle.  Compared  with  the  parti¬ 
cle  center  temperature  shown  in  Fig.  3,  the  average  temperature 
demonstrates  almost  the  same  rising  trend,  except  for  the  short 
lag  (a  few  seconds)  at  the  beginning.  With  the  furnace  wall  tem¬ 
perature  rising,  the  maximum  particle  temperature  rise  rate 
increased  from  9.7  Ks'1  (1000  K)  to  26.9  Ks"1  (1200  K)  due  to 
the  enhancement  of  the  radiative  heating  from  the  furnace.  How¬ 
ever,  the  time  corresponding  to  the  maximum  temperature  rise 
rate  remained  at  3  s  for  each  condition,  as  it  is  determined  by  both 
the  radiative  heating  and  intra-particle  heat  conduction. 

4.2.  Pyrolysis  of  straw  particles 

Figs.  7  and  8  illustrate  the  time  history  of  the  particle  center 
temperature  and  residual  mass  of  straw  particles  under  different 
furnace  wall  temperatures  (1000  K,  1100K,  and  1200  I<),  respec¬ 
tively.  All  particles  had  the  similar  initial  volume  (diameter  of 
8  mm)  and  initial  mass  (~360  mg).  Similar  to  the  coal  particle 
experiments,  each  experiment  was  repeated  three  times  and  the 
average  result  is  shown  together  with  the  error  bar  indicating 
the  statistical  uncertainty  of  the  three  measurements  at  each  con¬ 
dition.  As  illustrated  by  Fig.  7,  the  experimental  and  modeled  par¬ 
ticle  center  temperatures  of  the  straw  particle  show  a  considerable 
disagreement,  especially  at  the  early  stage  of  the  pyrolysis.  This 
disagreement  is  probably  due  to  the  swelling  of  straw  particles 
during  the  pyrolysis  process.  According  to  the  observation  during 
experiments,  the  volume  of  the  straw  particle  rapidly  expanded 
to  more  than  twice  larger  than  its  original  volume  and  became  a 
porous  structure  during  the  pyrolysis.  In  this  case,  the  thermocou¬ 
ple  may  not  have  maintained  contact  with  the  particle  interior. 
Owing  to  the  effect  of  radiative  heating  from  the  furnace,  temper¬ 
ature  measured  by  the  thermocouple  was  higher  than  the  real 
value  to  a  considerable  extent.  When  the  pyrolysis  finished,  the 
swelling  of  the  straw  particle  stopped  and  the  measured  tempera¬ 
ture  became  accurate  again. 

As  shown  in  Fig.  8,  the  bio-CPD  model  can  predict  the  pyrolysis 
of  straw  particles  under  different  temperatures  reasonably.  At  the 
time  of  200  s,  the  measured  total  volatile  yields  of  straw  particles 
were  71.4%,  72.2%,  and  71.6%  (dry  basis)  with  the  furnace  wall  tem¬ 
perature  at  1000  K,  1100  K,  and  1200  K,  respectively.  Since  the  vol¬ 
atile  mass  fraction  of  the  straw  in  the  proximate  analysis  was 
72.04%  (dry  basis;  measured  at  1173  K;  see  Table  1),  it  can  be  con¬ 
cluded  that  the  straw  particle  already  reached  the  maximum  yields 
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Fig.  8.  Time  history  of  residual  mass  of  straw  particles  under  different  furnace  wall 
temperatures. 


at  1000  K,  and  the  variance  of  the  yields  under  different  tempera¬ 
tures  was  due  to  experimental  error.  The  model  prediction  of  the 
total  volatile  yields  were  67.6%,  68.4%,  and  68.5%  (dry  basis)  with 
the  furnace  wall  temperature  at  1000  K,  1100K,  and  1200  K, 
respectively,  which  matched  the  measured  data  within  3.8  wt% 
upon  200  s  of  the  pyrolysis  process. 

The  time  history  of  the  mass  loss  rate  of  straw  particles  under 
different  furnace  wall  temperatures  is  shown  in  Fig.  9.  From  the 
experimental  data,  the  maximum  pyrolysis  rate  of  all  the  three 
cases  was  around  1.9%  s-1,  which  was  much  larger  than  that  of 
the  coal.  The  reason  lies  in  the  dynamic  equilibrium  of  radiative 
heating  and  the  cooling  of  pyrolysis.  While  the  furnace  tempera¬ 
ture  rose,  the  straw  particle  tended  to  devolatilize  faster;  however, 
faster  pyrolysis  could  enhance  the  cooling  effect  of  pyrolysis,  since 
the  pyrolysis  process  was  endothermic  and  the  yielded  gas  could 
help  cooling  the  particle.  As  also  shown  in  this  figure,  the  time  span 
of  the  maximum  pyrolysis  rate  was  relatively  long  (more  than 
25  s).  When  the  temperature  rose  from  1000  I<  to  1200  K,  the  time 
to  reach  the  maximum  pyrolysis  rate  decreased  from  about  66  s  to 
43  s.  Compared  with  the  experimental  data,  although  the  mass  loss 
rate  of  the  case  of  1200  K  was  overestimated,  the  bio-CPD  model 
predicted  the  trends  of  the  mass  loss  rate  at  an  acceptable  accu¬ 
racy.  It  should  be  stressed  that  there  is  no  adjustable  parameter 
in  the  model. 
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Fig.  7.  Time  history  of  center  temperature  of  straw  particles  under  different  furnace 
wall  temperatures. 


Fig.  9.  Time  history  of  mass  loss  rate  of  straw  particles  under  different  furnace  wall 
temperatures. 
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The  mass-weighted  average  temperature  of  straw  particles  and 
the  rise  rate  of  the  temperature  calculated  by  the  bio-CPD  model 
are  shown  in  Fig.  10.  Similar  to  the  results  of  coal  particles,  the 
mass-weighted  average  temperature  of  straw  particles  also  dem¬ 
onstrates  the  same  trend  (except  for  the  short  lag  at  the  begin¬ 
ning),  compared  to  the  particle  center  temperature  (see  Fig.  7). 
With  the  furnace  wall  temperature  rising,  the  maximum  particle 
temperature  rise  rate  increased  from  9.0 1<  s  1  (1000  K)  to 
22.1  K  s_1  (1200  I<),  which  was  a  bit  lower  than  that  of  the  coal  par¬ 
ticle.  Since  the  straw  yielded  more  volatile  than  the  coal,  its  cooling 
effect  was  also  larger.  The  time  corresponding  to  the  maximum 
temperature  rise  rate  was  3  s  for  the  straw  particle  in  all  the  three 
cases,  which  is  the  same  as  that  in  coal  particle  pyrolysis.  Com¬ 
pared  with  that  of  the  coal  particle,  the  overall  temperature  rise 
rate  of  the  straw  was  faster  and  it  experienced  a  second  peak  after 
the  mass  loss  rate  reached  the  maximum.  It  should  be  noted  that 
the  heat  capacity  of  the  straw  particle,  which  had  a  great  influence 
on  the  temperature  rise  rate,  decreased  faster  than  that  of  the  coal 
particle,  since  its  pyrolysis  rate  and  mass  loss  rate  were  larger.  In 
the  early  stage,  the  temperature  rise  rate  of  the  straw  particle 
declined  as  the  mass  loss  rate  increased,  because  the  radiation  heat 
flux  of  the  particle  decreased  and  the  cooling  effect  of  pyrolysis 
gradually  enhanced.  Flowever,  when  the  mass  of  the  straw  particle 
declined  to  a  critical  value,  its  heat  capacity  became  sufficiently 
small  to  balance  the  decreased  radiation  heat  flux  and  pyrolysis 
cooling  effect,  the  temperature  rise  rate  stopped  further  decrease 
but  began  to  elevate.  At  the  final  stage,  since  the  particle  temper¬ 
ature  was  very  high,  the  temperature  rise  rate  decreased  again 
due  to  the  small  radiation  heat  flux. 


4.3.  Pyrolysis  of  blended  coal-straw  particles 


Figs.  11  and  12  illustrate  the  time  history  of  the  center  temper¬ 
ature  and  residual  mass  of  blended  coal-straw  particles  with  dif¬ 
ferent  mass  ratios  (coal/straw  =  80:20,  50:50,  20:80), 

respectively.  All  particles  had  the  similar  initial  volume  (diameter 
of  8  mm)  and  initial  mass  (—360  mg),  and  the  furnace  wall  temper¬ 
ature  was  1100  I<  for  all  the  three  cases.  Open  symbols,  solid  sym¬ 
bols,  and  lines  in  Figs.  11  and  12  are  experimental  data,  weighted 
average  of  foregoing  individual  coal  and  straw  experimental  data, 
and  model  predictions,  respectively.  Each  experiment  was  also 
repeated  three  times  and  the  average  result  is  shown  together  with 
the  error  bar  indicating  the  statistical  uncertainty  of  the  three  mea¬ 
surements  at  each  condition.  As  Fig.  1 1  shows,  the  experimental, 
weighted  average  and  modeled  center  temperature  of  the  blended 
particle  agreed  well  with  each  other. 
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Fig.  10.  Time  history  of  mass-weighted  average  temperature  (lines  with  symbols) 
and  temperature  rise  rate  (lines)  of  straw  particles  under  different  furnace  wall 
temperatures. 
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Fig.  11.  Time  history  of  center  temperature  of  blended  coal-straw  particles  at 
different  mixing  ratios. 
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Fig.  12.  Time  history  of  residual  mass  of  blended  coal-straw  particles  at  different 
mixing  ratios. 


As  illustrated  by  Fig.  12,  the  residual  mass  fraction  of  coal-straw 
blends  calculated  through  the  weighted  average  of  foregoing  indi¬ 
vidual  coal  and  straw  experimental  data  also  agrees  well  with  the 
experimental  data  of  blends  with  different  ratios.  The  RMS  values 
of  the  relative  errors  vary  between  0.014  and  0.051,  suggesting 
therefore  the  absence  of  synergistic  effects  between  the  two  fuels 
during  the  pyrolysis  process.  This  result  is  also  supported  by  previ¬ 
ous  studies  of  co-pyrolysis  of  coal-biomass  blends  [7-11].  At  the 
time  of  500  s,  the  measured  total  volatile  yields  of  blended  particles 
were  38.3%,  51 .4%,  and  64.9%  (dry  basis)  with  the  coal/straw  ratio  of 
80:20,  50:50,  and  20:80,  respectively.  Compared  with  the  mass 
weighted  averages  of  the  volatile  mass  fraction  in  the  proximate 
analysis  of  the  coal  and  the  straw  (see  Table  1),  which  were  39.1%, 
51.5%,  and  63.8%  (dry  basis)  for  the  blended  particle  at  the  three 
mixing  ratios,  the  difference  is  less  than  1.1  wt%.  This  perfect  linear 
relationship  between  the  coal/straw  ratio  and  the  final  volatile 
yields  corroborates  that  there  is  no  interaction  between  the  coal 
and  the  straw  during  the  pyrolysis  process.  The  predicted  final  vol¬ 
atile  yields  using  the  extended  CPD  model  were  39.1%,  50.4%,  and 
61.4%  (dry  basis)  for  the  blended  particle  at  the  three  mixing  ratios, 
which  matched  the  measured  data  within  3.5  wt%.  The  model  pre¬ 
diction  fitted  the  experimental  results  in  all  different  cases,  which 
also  implies  the  effect  of  co-pyrolysis  ratio  is  insignificant.  The  fact 
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Fig.  13.  Time  history  of  mass  loss  rate  of  blended  coal-straw  particles  at  different 
mixing  ratios. 


5.  Conclusions 

Pyrolysis  of  coal,  straw,  and  blended  coal-straw  particles  were 
studied  in  a  single-particle  reactor  system,  with  the  temperature 
and  coal/straw  mixing  ratio  as  the  main  parameters.  The  statistics 
of  pyrolysis  of  the  coal  and  straw  agreed  well  with  the  experimen¬ 
tal  data  of  the  blends  at  different  coal/straw  mass  ratios,  indicating 
a  lack  of  synergistic  effects. 

A  one-dimensional,  time-dependent  single  particle  pyrolysis 
model  was  developed  based  on  the  chemical  percolation  devolatil¬ 
ization  (CPD)  model  and  the  upgraded  bio-CPD  to  simulate  the 
pyrolysis  of  coal  and  straw  particles.  Model  predictions  at  three 
pyrolysis  temperatures  agreed  well  with  the  experimental  data. 

An  extended  CPD  model  was  proposed  to  describe  the  co-pyro- 
lysis  characteristics  of  coal-biomass  blends.  Pyrolysis  of  the  blends 
at  different  coal/straw  mass  ratios  was  predicted  properly  by  the 
proposed  model  without  adjusting  any  model  parameters.  Encour¬ 
aged  by  the  agreement  between  the  model  prediction  and  the 
experimental  data,  further  work  will  be  done  to  expand  the  appli¬ 
cability  of  the  extended  CPD  model  to  other  coal-biomass  blends. 
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Fig.  14.  Time  history  of  mass-weighted  average  temperature  (lines  with  symbols) 
and  temperature  rise  rate  (lines)  histories  of  blended  coal-straw  particles  at 
different  mixing  ratios. 


that  the  macromolecular  structures  of  the  coal  and  the  straw  are 
cognate  carbon-molecular  frameworks  can  be  a  main  reason  for 
the  lack  of  interactions  between  the  two  fuels.  The  active  radicals 
produced  from  both  fuels  are  similar  and  there  is  no  catalytic  agent 
to  interfere  with  the  chemical  reactions  between  them,  under  the 
inert  atmosphere  [9  .  Therefore,  chemical  interaction  between  the 
two  fuels  is  hardly  observed. 

The  time  history  of  the  mass  loss  rate  of  blended  particles  at  dif¬ 
ferent  mixing  ratios  is  shown  in  Fig.  13.  The  measured  maximum 
pyrolysis  rates  rose  from  0.62%  s_1  to  1.91%  s1,  as  the  mass  frac¬ 
tion  of  straw  in  the  blend  increased  from  20%  to  80%.  However, 
the  time  corresponding  to  the  maximum  pyrolysis  rate  remained 
at  about  46  s  for  the  three  cases,  since  it  was  mainly  determined 
by  the  heating  condition.  In  comparison  with  the  experimental 
data,  the  extended  CPD  model  predicted  the  trends  of  the  mass  loss 
rate  reasonably. 

The  mass-weighted  average  temperature  of  blended  particles 
and  rise  rate  of  the  temperature  calculated  by  the  extended  CPD 
model  are  shown  in  Fig.  14.  Since  the  heating  condition  of  the  three 
cases  was  the  same,  the  curves  of  the  average  particle  temperature 
and  its  rise  rate  overlapped  in  most  of  the  pyrolysis  period.  The 
average  particle  temperature  tended  to  rise  quicker  as  the  mass 
fraction  of  straw  in  the  blend  increased. 
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